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ABSTRACT 

The hydrodynamical modelling of Type II plateau supernova light curves predicts a 
correlation between three observable parameters (the plateau duration, the absolute 
magnitude and photospheric velocity at the middle of the plateau) on the one side 
and three physical parameters (the explosion energy E, the mass of the envelope 
expelled M, and the presupcrnova radius R) on the other side. The correlation is 
used, together with adopted EPM distances, to estimate E, A4, and R for a dozen 
of well-observed SNe IIP. For this set of supernovae, the resulting value of E varies 
within a factor of 6 (0.5 < i?/10 51 erg < 3), whereas the envelope mass remains within 
the limits 10 < M/Mq < 30. The presupernova radius is typically (200 — 600) Rq, 
but can reach > 1000 Rq for the brightest supernovae (e.g., SN 1992am). 

A new method of determining the distance of SNe IIP is proposed. It is based 
on the assumption of a correlation between the explosion energy E and the 56 Ni 
mass required to power the post-plateau light curve tail through 56 Co decay. The 
method is useful for SNe IIP with well-observed bolometric light curves both during the 
plateau and radioactive tail phases. The resulting distances and future improvements 
are discussed. 

Key words: supernovae: general - galaxies: distances and redshifts. 



X 



1 INTRODUCTION 

Plateau Type II supernovae (SNe IIP) are believed to come 
from the explosion of massive supergiant stars whose en- 
velopes are rich in hydrogen. Their light curves are easy 
to identify by a long plateau (sometimes up to 120-150 
d) which is the result of the propagation of a cooling-and- 
recombination wave (CRW) through the supernova envelope 
that is in a state of free inertial expansion (u = r/t). The 
CRW physics is discussed in detail by Imshennik & Nady- 
ozhin (1964), Grassberg, Imshennik & Nadyozhin (1971), 
and Grassberg & Nadyozhin (1976). The CRW propagates 
supersonically downward through the expanding supernova 
envelope and separates almost recombined outer layers from 
still strongly ionized inner ones. During the plateau phase, 
the photosphere sits on the upper edge of the CRW front. 
Since the CRW downward speed turns out to be close to the 
velocity of the outward expansion, the photospheric radius 
changes only slowly during the plateau phase. If one takes 
into account that also the effective temperature does not 
change appreciably (it approximately equals the recombina- 
tion temperature 5000-7000 K), the approximate constancy 
of the luminosity becomes obvious. 



The supernova outburst properties are determined 
mainly by three physical parameters: the explosion energy 
E, the mass M of the envelope expelled, and the initial ra- 
dius R of the star just before the explosion (presupernova). 
Litvinova & Nadyozhin (1983, 1985) have undertaken an at- 
tempt to derive these parameters from a comparison of the 
hydrodynamical supernova models with observations. They 
constructed simple approximation formulae which allow to 
estimate E, M, and R from the observations of individual 
SNe IIP. Their results were confirmed by an independent 
semi-analytical study (Popov 1993). At that time, only one 
or two supernovae were sufficiently observed to apply these 
formulae. At present, there exist detailed observational data 
for 14 such supernovae, including in 12 cases expanding pho- 
tosphere (EPM) distances, which we use in section 2 to es- 
timate E, M, and R by means of these formulae. 

In section 3, we propose a new method of distance de- 
termination and employ it to 9 individual SNe IIP which 
are well-observed both at the plateau and radioactive-tail 
phases. The method is based on the assumption of a corre- 
lation between the explosion energy E and the mass of 56 Ni 
in the supernova envelope. In section 4 we compare phys- 
ical parameters and distances of SNe IIP as derived from 
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Figure 1. A schematic SN IIP light curve. An open circle marks 
the middle of the plateau and two black circles show the plateau 
boundaries. The light curve tail powered by the 56 Co decay is 
also shown (tc = 111.3 d). 

the new method with those obtained previously from the 
EPM method and discuss also other aspects of our results. 
Concluding remarks are given in section 5. 

The preliminary results of this study were reported to 
the Workshop on physics of supernovae held at Garching, 
Germany, July 2002 (Nadyozhin 2003). 



2 A COMPARISON OF HYDRODYNAMIC 
MODELS WITH OBSERVATIONS 

Figure 1 shows a schematic SNellP light curve. The plateau 
is defined as part of the light curve on which the super- 
nova brightness remains within 1 mag of the mean value. For 
some supernovae, the plateau begins almost immediately af- 
ter the onset of the explosion (t — 0) whereas for others a 
short luminosity peak can precede the plateau. The peak ei- 
ther appears owing to a shock wave breakout in the case of 
presupernovae of not very large initial radii (R < 1000 Rq) 
or, according to Grassberg et al. (1971), originates from the 
emergence of a thermal wave precursor for presupernovae of 
very large radii (R « (2000 — 5000) Rq) and of moderate 
explosion energies (E < 1 x 10 51 erg), or at last it may occur 
as a result of interaction between the supernova envelope 
and a dense stellar wind (Grassberg & Nadyozhin 1987). 
For some SNe IIP the peak duration 5t lasts only a few days 
and is difficult to observe (shock wave breakout), for others 
it could be as large as 10-20 days (thermal wave or dense 
wind) - examples for the latter may be such supernovae as 
SNel988A, 1991al, and 1992af (see below). 

It is quite clear that the middle of the plateau is to 
be used as the main reference point to compare the the- 
oretical models with observations. Litvinova & Nadyozhin 
(1983, 1985) (LN83 and LN85, hereafter) calculated a grid 
of the supernova models for E, Ai, and R within limits of 
(0.18 - 2.91) x 10 51 erg, (1 - 16) Mq, and (300 - 5000) Rq. 
They found E, M, R to be strongly correlated with the 
plateau duration At, and the middle-plateau values of the 
absolute ^-magnitude My, and the expansion velocity u p h 
at the level of the photosphere (Fig. 1). According to LN85, 
the following approximate relations can be used to derive E, 
M, and R from observations: 

\gE = 0.135M v + 2.341gAt + 3.131gitp h - 4.205, (1) 



IgM = 0.234 M v + 2.91 lg At + 1.96 lgw ph - 1.829, (2) 
\gR = -0.572 M v - 1.07 lg At- 2.74 lgw ph - 3.350, (3) 

where E is expressed in units of 10 51 erg, M. and R are in 
solar units, At in days, and u p h in 1000 km s -1 . Here Mv 
can be expressed through the apparent V magnitude by the 
relation: 

M v = V -A v - 51g(D/lMpc) - 25, (4) 

where D is the distance to a supernova and Av is the total 
absorption on the way to the supernova. One can find from 
Eqs. (l)-(3) that E, M, and R scale with the distance as: 

E~D- 0S75 , M ~ D~ 117 , R ~ D 2 - 86 . ( 5 ) 

Thus, it is very important to know D with as high accuracy 
as possible. We have selected 14 SNe whose observational 
data are collected in Table 1. The entries are: the heliocen- 
tric recession velocities v (from the NED: NASA/IPAC Ex- 
tragalactic Database) in column 3, the total absorption Av 
in column 4, the apparent V magnitude of the mid-point 
of the plateau in column 5, the duration At of the plateau 
in column 6, and the photosphere expansion velocity u p h in 
column 7. The references are in column 8. 

In order to check the extrapolative capability of 
Eqs. (l)-(3), we have included the SN 1987A in our anal- 
ysis. It is well known, that the SN 1987A presupernova ra- 
dius was as small as « 50 Rq - i.e. outside the interval of 
(300 — 5000) Rq encompassed by the above equations. More- 
over, the major part of the SN 1987A plateau (about 70 of 
110 d) was powered by the 56 Co-decay (see the review of 
Imshennik & Nadyozhin 1989 and references therein). 

Derived properties of the 14 SNe IIP are in Table 2. 
Column 2 is the recession velocity V220 of the supernova 
(column 1) corrected for a self-consistent Virgocentric in- 
fall model with a local infall vector of 220 km s -1 as de- 
scribed by Kraan-Korteweg (1986). Column 3 gives the 
distance Dh = V220/H0 assuming arbitrarily a value of 
Ho — 60kms _1 Mpc -1 . For comparison, column 4 gives the 
distance Depm obtained with the use of the expanding pho- 
tosphere method in the references listed at the bottom of the 
table. The SNe 1991al and 1992af are the exception. Owing 
to the incompleteness of the observational data, it is hard 
to determine the EPM distance to the SN 1991al (Hamuy 
2001). For the same reason, the EPM distance of 55 Mpc 
for the SN1992af obtained by Schmidt et al. (1994a) seems 
to be quite uncertain as pointed out by Hamuy (2001). For 
these two SNe, we present in column 4 the distances cal- 
culated by Hamuy (2001) from the CMB redshifts and the 
Hubble constant Ho = 65kms~ 1 Mpc -1 . Columns 5-8 are 
the absolute magnitude Mv of the mid-point of the plateau, 
the explosion energy E, the mass expelled Ai, and the presu- 
pernova radius R - all derived from Eqs. (4), (l)-(3) for the 
Dh distances listed in column 3. Column 9 gives the mass 
of 56 Ni, ejected by some supernovae, which was estimated 
by reducing the radioactive-tail luminosities, measured by 
Hamuy (2001), to the distances Dh given in column 3. 

For SN1987A, the resulting values of E and M 
(Table 2) differ no more than by a factor of 1.5 from cur- 
rent estimates based on a detailed study. However, the pre- 
supernova radius turned out to be too large. This happened 
because the LN85 approximations do not take into account 
the radioactive heating. An advanced study (Grassberg & 
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Table 1. Observational data for 14 SNcIIP. 



QTVT 


Host galaxy 


i>0 


Ay 


V 


at 


« P h 


Hei. 






km s — 1 


mag 


mag 


days 


km s — 1 




W 


( 2 ) 


(3) 


(4) 


(5) 


(6) 




(8) 


1968L 


NGC 5236 


516 


0.219 


12.0 


80 


4100 


1, 2, 3 


1969L 


NGC 1058 


518 


0.203 


13.4 


100 


4000 


1, 2 


1986L 


NGC 1559 


1292 


0.099 


14.7 


110 


4000 


4 


1988A 


NGC 4579 


1519 


0.136 


15.0 


110 


3000 


1, 2, 4, 5, 6 


1989L 


NGC 7339 


1313 


1.00 


16.5 


140 


3000 


7, 19 


1990E 


NGC 1035 


1241 


1.083 


16.0 


120 


4000 


2, 4, 8, 9 


1991al 


LEDA 140858 


4572 


0.318 


17.0 


90 


6000 


4 


1992af 


ESO 340-G038 


6000 


0.171 


17.3 


90 


6000 


4, 7 


1992am 


anon 0122-04 


14600 


0.464 


19.0 


110 


4800 


4, 10 


1992ba 


NGC 2082 


1104 


0.193 


15.43 


100 


2900 


4, 7 


1999cr 


ESO 576-G034 


6069 


0.324 


18.6 


100 


3600 


4 


1999em 


NGC 1637 


717 


0.314 


14.0 


110 


3000 


4, 11, 12, 13, 17 


1999gi 


NGC 3184 


592 


0.65 


15.0 


110 


2900 


14, 15, 16, 18 


1987A 


LMC 


278 


0.465 


3.3 


110 


2900 


4 



References. - (1) Patat et al. 1993; (2) Schmidt ct al. 1992; (3) Wood & Andrews 1974; 

(4) Hamuy 2001; (5) Ruiz-Lapuente et al. 1990; (6) Turatto et al. 1993; (7) Schmidt et al. 

1994a; (8) Schmidt et al. 1993; (9) Benetti et al. 1994; (10) Schmidt et al. 1994b; 

(11) Hamuy ct al. 2001; (12) Haynes et al. 1998; (13) Baron et al. 2000; (14) Schlcgcl 2001; 

(15) Smartt ct al. 2001; (16) Li ct al. 2002; (17) Elmhamdi ct al. 2003; 

(18) Leonard et al. 2002b; (19) Pennypacker & Pcrlmuttcr 1989 



Table 2. The supernova physical properties. 



SN 


t>220 


D H 


Depm 


M v 


E 


M 


R 


Mmo 




kms -1 


Mpc 


Mpc 


mag 


10 jl erg 


Mq 


Rq 


Mq 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


1968L 


291 


4.85 


4.5« 


-16.65 


0.83 


10.3 


286 




1969L 


766 


12.77 


10.6W 


-17.33 


1.05 


13.0 


595 




1986L 


1121 


18.68 


16.0M 


-16.76 


1.56 


23.5 


251 


0.026 


1988A 


1179 


19.65 


20.0W 


-16.60 


0.67 


14.5 


452 


0.082 


1989L 


1556 


25.93 


17.0« 


-16.57 


1.18 


29.8 


334 




1990E 


1238 


20.63 


18.0M 


-16.66 


1.98 


31.9 


200 


0.052 


1991al 


4476 


74.60 


70.0< 2a ) 


-17.68 


2.61 


17.6 


347 


0.12 


1992af 


6000 


100.00 


83.70< 2a ) 


-17.87 


2.46 


15.9 


445 


0.24 


1992am 


14600 


243.33 


180.0' 1 ) 


-18.40 


1.66 


13.9 


1321 


0.36 


1992ba 


1096 


18.27 


22.0^ 


-16.07 


0.57 


13.7 


272 


0.029 


1999cr 


6069 


101.15 


86.0( 2 ) 


-16.75 


0.90 


14.5 


368 


0.085 


1999em 


743 


12.38 


8.20< 3 ) 


-16.78 


0.63 


13.2 


569 


0.058 


1999gi 


707 


11.78 


11.10< 4 ) 


-16.01 


0.72 


18.7 


226 


0.025* 


1987A 




0.05 


0.05 


-15.66 


0.80 


22.6 


143 


0.065 



(^Eastman, Schmidt, & Kirshncr 1996; ( 2 >Hamuy 2001; ( 2 °'Hamuy 2001, based 

on an adopted value of Ho = 65, see text; ( 3 ) Leonard et al. 2002a; 
< 4 ) Leonard et al. 2002b 

*Derived from V = 17.86 mag at t = 174.3 d (Leonard et al. 2002b), and Hamuy's 



recipe (section 5.3 of his Thesis) to convert V into luminosity L. 



Nadyozhin 1986) demonstrates that the radioactive heating 
influences only weakly on E and M. furnished by Eqs. (1)- 
(3), whereas the R- values can be overestimated by a factor 
of 3. In this connection, one should have in mind that for 
some supernovae the R- values from Table 2 can be larger 
than actual presupernova radii. 

According to Table 2, the resulting values of E, M, 
and R seem to be reasonable enough: the expelled mass, 



explosion energy, and presupernova radius remain approx- 
imately in limits (1O-3O)7W , (0.6 - 2.6) x 10 51 erg, and 
(200 — 1300) -Rq, respectively. Hamuy (2001) assumed that 
SNe 1991al, 1992af were discovered several weeks after the 
explosion. Their plateaus, therefore, could have lasted for 
At ~ 110 d. It is quite probable, however, that their peak 
duration was St w 20 d for the reasons mentioned above. 
Having this in mind, we have chosen in Table 1 At = 90 d 
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which results in u p h = 6000 km s _1 . In the case of At = 110 d 
we would have to assume u p h = 7000 km s -1 and would ob- 
tain very large values of E and M for both supernovae: 
E « 7xl0 51 erg and M ~ 40A1q. No other special ad- 
justments of the observational data given in Table 1 were 
made. 



3 PLATEAU- TAIL DISTANCE 
DETERMINATION 

The SN IIP light curve tails are believed to be powered by 
the 56 Co decay. The temporal behavior of the bolometric 
luminosity is given by (see, e.g. Nadyozhin 1994): 

L = 1.45X10 43 exp (-—) erg s" 1 , (6) 

V TCo/ M.Q 

where t is measured from the moment of explosion (t — 0), 
AlNio is the total mass of 56 Ni at t = which decays with a 
half-life of 6.10 d into 56 Co, and t Co = 111.3 d. 
Equation (6) can be written in the form 

Mmo = ^Q, Q = F 41 (t)exp( TI ^), (7) 

where A^Nio is in Mq, t in days and D in Mpc. The quantity 
F41 (t) is the bolometric tail luminosity measured at time t in 
units 10 41 ergs -1 under the assumption that the supernova 
is at distance D — 1 Mpc. Equation (7) contains a single 
observational parameter Q which is independent of t and 
also makes no assumption on D as long as F41 is fixed by 
observations. 

Thus, it is irrelevant at which t the luminosity is actu- 
ally measured - one has only to be sure that the supernova 
really entered its tail phase. Columns 8 — 10 of Table 3 give 
t and corresponding values of F±i(t) and Q derived from 
Hamuy's (2001) Figures 5.7 and 5.8 except SN 1999gi for 
which the values were calculated from the data of Leonard 
et al. (2002b). 

If the value of .Mi-ho was known, one could easily find 
the distance D from Eq. (7) . So, we have to look for a way to 
estimate Mmo independently. It seems reasonable to assume 
that the supernova explosion energy E should correlate with 
Almo produced during the explosion. This means that 

E = f(M KI0 ) = f(j^Cj), (8) 

where / represents a statistically admissible correlation func- 
tion rather than a strict mathematical relation. Inserting 
this expression for E into Eq. (1) and using Eq. (4) for My, 
we obtain an equation which can be solved for D when 
V — Ay, w p h, At, and Q are known from observations. Then 
for given D, we can find E, M, R, and A^Nio from Eqs. (1)- 
(3), Eq. (4), and Eq. (7), respectively. 

What can be said about the function / (Mmo) at 
present, when the details of the SNII mechanism remain 
still ambiguous? First of all, it is reasonable to assume that 
a good fraction of E comes from the recombination of free 
neutrons and protons into 56 Ni just at the bottom of the en- 
velope to be finally expelled (Nadyozhin 1978, Bethe 1996). 
The hydrodynamical modelling of the collapse (Nadyozhin 
1978) have indicated that under favourable conditions a 
neutron-proton shell could be accumulated just under the 



steady accreting shock wave. When the mass of such a 
shell reaches some critical value (presumably of the order of 
wO.IA^q) the shell can become unstable in respect to re- 
combining into the "iron group" elements (specifically into 
56 Ni) to supply the stalled shock wave with the energy of 
?s 10 5 erg necessary to trigger the supernova. Here, there 
is a physical analogy with the origin of planetary nebulae 
from red giants where the energy from the recombination 
of hydrogen and helium causes the expulsion of a red giant 
rarefied envelope. The recent study (Imshennik 2002, and 
references therein) of the "neutrino crown" - the region en- 
closed within neutrinosphere and accreting shock, turns out 
to be in line with such a picture of the supernova mechanism 
. However, some Ni can be produced through the explosive 
carbon-oxygen burning induced by the outgoing shock wave. 
In this case the energy release per unit Ni mass is lower by 
an order of magnitude than for the neutron-proton recom- 
bination. 

The energy released by the neutron-proton recombina- 
tion, producing a 56 Ni mass of Mmo, is given by 

S(np-^Ni) = 1.66 xlO 52 erg. (9) 

Mq 

Thus, the production of only ~ 0.06 Mq of 56 Ni is suffi- 
cient to provide the standard explosion energy of 10 51 erg. 
The current hydrodynamic models of the SNII explosions 
(Woosley & Weaver 1995; Rauscher et al. 2002) do not show 
a correlation between E and A^Nio because in these models 
56 Ni comes from explosive silicon and carbon-oxygen burn- 
ing near to the envelope bottom and its yield is sensitive 
to the mass cut point. The photometrical and spectroscop- 
ical properties of the SN models are virtually independent 
of the mass cut. On the contrary, the nucleosynthesis yields 
are very sensitive to the mass cut. In the current SN models 
the explosion is usually simulated by locating a piston at the 
internal boundary m = Aleut- The piston moves with time 
according to a prescribed law, 7? p j s (t), with the velocity (R) 
amplitude being chosen to ensure the final kinetic energy of 
the expelled envelope of the order of 10 jl erg. There are two 
major uncertainties at this point. First, for a given velocity 
amplitude the resulting nuclear yields are still sensitive to 
the form of the function R p i B (t). Second, the presupernova 
structure (especially chemical composition) in the vicinity 
of m = Mcut will always remain ambiguous until the de- 
tailed mechanism of the SN disintegration onto the collapsed 
core and thrown envelope is established. The point is that 
such 2D effects as rotation and large-scale mixing can result 
in the presupernova structure different from that predicted 
by the spherically symmetrical models. Under such circum- 
stances, it is difficult to find a serious argument against the 
possibility to expel a noticeable amount of 56 Ni from the re- 
combination of the neutron-proton shell. Thus, we propose a 
neutron-proton layer which is located somewhat deeper than 
the value of Af cu t assumed in the current SN models. This 
layer recombines into 56 Ni providing the energy sufficient 
to convert a steady-state accretion shock into the outgoing 
blast wave. In this case a good correlation between E and 
AIniO is to be expected. 

The proposed correlation can have a complex nature. 
It is quite probable that the function / in Eq. (8) depends 
also on M since the supernova mechanism is expected to be 
sensitive to the presupernova mass. For us only the existence 
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Table 3. The tail-calibrated supernova physical properties (£ = 1). 



SN 


-Dp-T 


M v 


E 


M 


R 


Mmo 


F 4 i(t) 


t 


Q 




Mpc 


mag 


10 51 erg 


Mq 


Rq 


Mq 


10 41 erg/(s Mpc 2 ) 


days 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


1986L 


29.67 


-17.76 


1.14 


13.7 


944 


0.067 


2.25 xl0~ 3 


180 


0.0113 


1988A 


15.21 


-16.05 


0.79 


19.6 


217 


0.048 


4.96 xl0~ 3 


200 


0.0299 


1990E 


29.16 


-17.41 


1.57 


21.3 


539 


0.094 


2.67 xl0~ 3 


200 


0.0161 


1991al 


85.31 


-17.97 


2.38 


15.0 


509 


0.14 


8.13xl0~ 4 


140 


0.00286 


1992af 


86.45 


-17.55 


2.71 


18.8 


293 


0.16 


9.02X10- 4 


140 


0.00317 


1992ba 


19.85 


-16.25 


0.53 


12.4 


346 


0.032 


1.97 xl0~ 3 


200 


0.0119 


1999em 


11.08 


-16.54 


0.68 


15.0 


414 


0.041 


1.26 xl0~ 2 


150 


0.0485 


1999gi 


14.53 


-16.46 


0.63 


14.5 


411 


0.038 


5.41 xl0~ 3 


174 


0.0259 


1987A 


0.045 


-15.42 


0.87 


25.6 


104 


0.053 


8.16X10 2 


170 


3762 



of some correlation is important which in combination with 
Eqs. (l)-(3) allows to determine the distance independently. 

To demonstrate how such a method can work we make 
the simplest assumption that E is proportional to D(np — > 
Ni). Then one can write: 

E = £D(np^Ni) = 16.6£A4 N io = 0.1145 £D 2 Q , (10) 

where, as usual, E is in 10 51 erg, M Nio in A4q and D in 
Mpc. This equation implies that the function /, introduced 
in Eq. (8), reads as f(x) = 16.6 £ x where £ is an adjustable 
parameter which can be either less or larger than 1. If there 
is a noticeable contribution to A^jvio from the explosive 
carbon-oxygen burning then £ < 1; if a noticeable contribu- 
tion to the explosion energy comes from other source rather 
than the neutron- proton recombination then £ > 1. 

Inserting E from Eq. (10) and My from Eq. (4) into 
Eq. (1) and solving for D, we obtain: 

lgD = -0.374 lg(£Q) + 0.0504 (V - A v ) + 0.875 lg At 
+1.17 lgwph- 2.482, (11) 

where D is in Mpc, At in days, and u p h in 1000 km s^ 1 . We 
will refer to distances derived from Eq. (11) as 'plateau-tail 
distances', Dp-T, hereafter. The results are given in Table 3 
for nine supernovae selected from Table 2. We did not in- 
clude SNe 1992am and 1999cr in our analysis because their 
last available observations may not yet reflect the radioac- 
tive tail phase. Specifically, there are only two observations 
of SN 1992am at the post-plateau phase of the light curve. 
Since the observations are separated by a short time interval 
of 3 days, it is difficult to derive the inclination of the bolo- 
metric light curve with a required accuracy to be shure that 
SN 1992am is already in the radioactive-tail phase. More- 
over, one has to remember that in addition to the Co-decay 
the tail luminosity can also be contributed by the ejecta- 
wind interaction (see Chugai 1991 and references therein). 
SN 1992am is suspicious in this respect because its presu- 
pernova radius seems to be larger than 1000 Rq (Table 2). 
Hence, the A4Nio-values for these SNe in Table 2 could be 
actually upper limits. 

The different columns of Table 3 give the following 
quantities: (2) the distance Dp-t from Eq. (11) setting 
£ = 1; (3) the corresponding absolute V-magnitude of the 
mid-point of the plateau Mv; (4)-(7) the quantities E, M, 
R, and A^Nio as in Table 2, but now using the distance Dp-t 
as in column (1); the columns (8)-(10) are explained above. 



The values of E, M, R, and A^mo for the ^-values differ- 
ent from 1 can be found using the following scaling relations 
which result from Eqs. (5), (7), and (11): 

E ~ e 252 , M ~ £ ' 438 , R ~ r 1 ' 07 , A^Nio ~ r°- 748 - (12) 

For a fixed Q, the dependence of the distance Dp-t, defined 
by Eq. (11), on extinction Av proves to be very weak: an er- 
ror in Av of ±1 mag changes Dp-T by only ±12%. However, 
if the tail luminosity F41 is derived from the V measurements 
(just the case of Hamuy's J4i-values we use here) then the 
lg F41, and consequently lg Q, scale as 0.4 Av and lg Dp-t, 
derived from Eq. (11), actually varies with Av in a standard 
way, as —0.2 Ay. If the tail luminosity were derived from 
infra-red measurements then the resulting Dp_T distances 
would be largely independent of extinction. Note also rather 
weak dependence on £Q: Dp-T ~ (£<3)~°' 374 - For instance, 
the decrease in £Q by a factor of 2 results in an increase of 
Dp-t by 30% only. 

The random errors typically of ±10% for the St and M p h 
values assumed in Table 1 result in the uncertainty factor 
of w 1.2 for Dp-t and « 1.5 for A4nh)(~ D 2 ) given in Ta- 
ble 3. However, one has to keep in mind two main sources 
of systematic errors: (i) probable deviation of the theoret- 
ical models (which Eqs. 1-3 are based on) from real SNe, 
and (ii) the presentation of the E — A^nm correlation in the 
form of the straight proportionality (Eq. 10). Both the types 
of systematic errors are difficult to estimate at present. Al- 
though the SN models calculated in LN83 and LN85 rest 
upon a very simplified presupernova structure, they consis- 
tently take into account the ionization and recombination of 
hydrogen and helium thereby remaining still useful. When 
a new grid of the SN models, based on modern evolution- 
ary presupernova structure, is created the systematic error 
(i) certainly will be reduced. The reduction of the system- 
atic error (ii) requires a more profound knowledge of the 
SN mechanism. Empirically, this problem can be solved by 
adjusting the factor £ for each individual SN. It is necessary, 
however, to collect a much more rich statistics (at least by 
a factor of 3) than that available nowadays (only 9 SNe in 
Table 3). 



4 DISCUSSION 

The plateau-tail distances derived in section 3 and listed in 
column 2 of Table 3 are plotted in a Hubble diagram in Fig. 2 
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Figure 2. The Hubble diagram of 8 SNellP with D P _ T dis- 
tances from plateau and tail observations (black circles). Also 
shown are the eleven SNellP with known EPM distances (as- 
terisks). The respective Hubble lines are fitted to the data. The 
abscissa is the distance modulus (m — M) = 5 lg D + 25. 



(except SN 1987A which is not in the Hubble flow) . The eight 
SNellP define a Hubble line with H = 55±5 kms" 1 Mpc -1 . 
Also shown in Fig. 2 are the eleven SNe IIP for which EPM 
distances have been published (column 3 of Table 2). They 
define a Hubble line of H = 70 ±4kms _1 Mpc -1 , i.e. the 
EPM distances are smaller than the plateau-tail distances 
by 25% on average. 

At this point it is not possible to decide which of the two 
results is more nearly correct. Both methods, the plateau- 
tail distances and the EPM distances, depend on assump- 
tions which are difficult to verify. The EPM method faces 
the problem of the dilution factor in an expanding atmo- 
sphere and the definition of the photospheric radius which 
depends on the uncertainties connected with the opacity of 
an expanding medium. However, it may be noted that the 
EPM distance of SN 1987A agrees well with the generally 
adopted distance of LMC of 50 kpc (Eastman, Schmidt, & 
Kirshner 1996) and the EPM distance of SN 1968L is indis- 
tinguishable from the Cepheid distance of NGC 5236 (M83) 
(Thim et al. 2003). 

The main assumption which affects the plateau-tail dis- 
tances concerns the nature of the proposed E — AAnm cor- 
relation. For our simplified example of such a correlation, 
all the uncertainties turn out to be cumulated in the pro- 
portionality factor £ between the explosion energy E and 
the nickel mass A4Nio- In Table 3 we have adopted a plau- 
sible value of £ = 1, but it cannot be excluded that £ is as 
low as 0.5 or as high as 2. Since the Hubble constant scales 
as Ho ~ £° ,374 , an average value as high as £ = 1.9 would 
be needed to bring the plateau-tail distances in general ac- 
cord with the EPM distances. Such a high average value 
of £ is, however, not supported by SNe 1987A and 1999gi. 
If the Dp-t distance of SN 1987A from Table 3 is scaled 
to the canonical LMC distance of 50 kpc, £ becomes 0.75. 
And if the host galaxy NGC 3184 of SN 1999gi with a D P _ T 



O 

CD o 
II 

X 



1 - 



1990E 



1986L 

o 



o E(H = 60) 
* E(EPM) 







1 



E(P-T)° 

Figure 3. The explosion energies E(H = 60) and B(EPM) 
versus E(P — T) (see the text). 
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Figure 4. The expelled masses M(Hq = 60) and A / !(EPM) 
versus M(P — T) (see the text). 



distance of 14.53 Mpc is a member of the same group as 
NGC 3198 and NGC 3319, for which Freedman et al. (2001) 
give a mean Cepheid distance of 13.5 Mpc, £ becomes 1.2. 
Eventually additional SNellP with large distances, where 
the influence of peculiar motions are negligible, will better 
determine the scatter of the Hubble diagram and allow a 
meaningful determination of the actual range of £. 

We have considered three sets of the physical supernova 
parameters E, Ai, and R: (i) for the Hubble distances Dh 
with H = 60kms _1 Mpc -1 (Table2, column 3); (it) for 
the EPM distances Depm (Table 2, column 4); (m) for the 
plateau-tail calibrated distances Dp-t (Table 3, column 2). 

Although the above parameters derived from the EPM- 
distances are not presented in Table 2, the corresponding E 
and .M-values can be read out of Figs. 3 and 4 which com- 
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Figure 5. The explosion energy-envelope mass diagram for the 
case of EPM distances -Depm from column 4 of Table 2 (aster- 
isks; SNc 1991al and 1992af being excluded) and for the case 
of the plateau-tail distances _Dp_T from column 2 of Table 3 
(black circles). Some SNe are identified (see text). 



pare E and M for sets (i) and (ii) with those for set (Hi). For 
seven SNe E and M are rather insensitive to the adopted 
distances. However for SNe 1986L and 1990E, labelled in 
Figs. 3 and 4, the deviations from the (P-T)-values are rather 
large, especially in case of the envelope mass M. These SNe 
differ from others by having a long plateau of (110 — 120) d 
in combination with still a substantial expansion velocity of 
4000 km s -1 . As a result, their envelope masses M, derived 
from the distances defined by the Dh and Depm values, ex- 
ceed those for other SNe. Such a discrepancy for these two 
SNe is considerably weakened if £ « 2. Such a high value 
of £ implies that half of the explosion energy is supplied by 
a source different from the neutron-proton recombination. 
This may indicate that for massive SNe the envelope mass 
M (in addition to Mnio) is involved in the correlation given 
by Eq. (8). 

The random errors of E and M from our approximate 
Eqs. (l)-(3) are estimated to be about ±30% . Observa- 
tional errors especially in the expansion velocity u p h and 
the plateau duration At can modify E and M by another 
factor of 1.3. Thus it seems reasonable to assume a random 
uncertainty of a factor of ~ 1.5 for the individual values of 
E and M in Tables 2 and 3. The presupernova radii R are 
very sensitive to distance errors (cf. Eq. 5) and may carry 
random errors of a factor of 2. The radii of SNe with large 
nickel masses like SN 1991al, 1992af and perhaps 1992am 
may carry additional systematic errors because Eqs. (l)-(3) 
do not take into account the radioactive heating in a consis- 
tent way. 

The expelled masses M are plotted against the explo- 
sion energies E in Fig. 5 for two cases, i.e. based on EPM 
and plateau-tail distances. In case of the Dp_t distances the 
mean mass of the eight SNe IIP is 16 .Mq with an rms devia- 
tion of only 3.Mq. This narrow mass range is contrasted by 
a wide range of explosion energies of (0.5 — 2.7) xlO 51 erg. 



The conclusion that there is no correlation between the ex- 
pelled mass - which is only (1.4 — 2) Ai q smaller than the 
presupernova mass - and the explosion energy is somewhat 
weakened by the values of Ai and E based on the EPM dis- 
tances suggesting a marginal correlation between Ai and E 
which is mainly due to only two SNe: 1986L and 1990E. 

One can think of a number of parameters which may 
explain the wide range of explosion energies. It could be ro- 
tation and magnetic fields inherited by the collapsing stel- 
lar core. It could be also nonspherical jet-like perturbations 
of a random nature arising from the macroscopic neutrino- 
driven advection below the accretion shock. Such pertur- 
bations could launch the outgoing blast wave earlier when 
the recombination nuclear energy stored in a hot neutron- 
proton gas was not yet as large as it should be in the case 
of spherical symmetry. If this is correct, one may expect 
that the asphericity of the explosion anticorrelates with the 
explosion energy. 

Recently, a promising project has been started (Van 
Dyk et al. 1999; Smartt et al. 2001, 2002; and references 
therein) with the ultimate aim to identify the supernova 
progenitors (presupernovae) or at least to impose conclu- 
sive constraints on their masses by inspecting the predis- 
covery field of nearby supernovae. In particular, Smartt et 
al. derived upper mass limits of 12 .Mq and 9A4q for the 
progenitors of the SNe 1999em and 1999gi, assuming dis- 
tances D for the host galaxies NGC 1637 and NGC3184 
of 7.5 Mpc and 7.9 Mpc, respectively. Note that these up- 
per limits depend on D and have to be adjusted for other 
values of D to 12 Mq(D/7.5 Mpc) ' 6 for SN 1999em and 
9.M Q (.D/7.9Mpc) ' 6 for SN 1999gi. This follows from the 
fact that the mass-luminosity relation can be approximated 
as L ~ A^ 3 ' 3 in the mass interval (10 — 15)^0. For 
SN 1999em at D P _ T = 11.08 Mpc (Table 3) follows 15.2 Mq 
as the upper mass limit for the SN 1999em progenitor. 
Hence, our result Ai = 15.0 Mq (Table 3) does not con- 
tradict to the observations as long as D(1999em) > 10 Mpc. 
The situation of SN 1999gi is similar. The upper mass 
limit for D(1999gi) = 14.53 Mpc (Table 3) is M < 9 x 
(14.53/7.9) 6 = 13.0 Mq, i.e. not in significant contradic- 
tion with the .M-value of 14.5 Mq from Table 3. There is no 
contradiction either with the upper mass limit of 15±| Mq 
for the SN 1999gi progenitor imposed recently by Leonard 
et al. (2002b). 

Equations (l)-(3) by LN85, derived from a grid of 
23 SNe IIP models covering a wide parameter space, im- 
ply a correlation between the absolute magnitude My (and 
hence luminosity L - both measured at the mid-point of the 
plateau) and the expansion velocity u p h. The correlation 
is shown in Figs. 6 and 7 where 23 grid models are shown 
by black dots; the straight lines are the least-squares fits. 
In Fig. 7 are also shown the eight observed SNe IIP from 
Table 3 marked by open circles, their absolute magnitudes 
My (Table 3, column 2) being calculated from Eq. 4, where 
the plateau-tail distances Dp_x were used from Table 3, 
column 2. These real SNe follow about the slope of the mod- 
els, but at a fixed value of tt p h they are fainter by « 0.6 mag 
on average. 

Empirically, Hamuy & Pinto (2002) have also found, 
using the CMB redshift-based distances, such a correlation. 
The slopes of their least-squares fits are virtually the same 
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for 23 SN models (black dots); open circles relate to 9 real SNc 
including SN 1987A. 



as shown in Figs. 6 and 7. Thus our models confirm their 
finding. 

The main conclusion one can draw from Figs. 6 and 7 
is that our three-parametric grid of only 23 SNellP prop- 
erly chosen models is ample enough to reproduce the main 
features of the real SNe. 



the plateau together with its total duration At are needed 
as additional input parameters. The results are presented in 
Table 2. 

Instead of using EPM or velocity distances it is also pos- 
sible to use the bolometric fluxes observed during the SN IIP 
tail phase to determine the Ni mass and hence new, indepen- 
dent distances called here plateau-tail distances Dp-t (cf. 
Eq. 11). The Dp-t distances yield new values of E, M, and 
R given in Table 3 for nine SNe which were observed both 
during their plateau and tail phases. The values of E and 
M, based on EPM and P-T distances agree well, with the 
exception of SNe 1986L and 1990E whose masses M com- 
ing from P-T distances are by a factor of 2 lower than from 
EPM distances (see Fig. 4) . 

The P-T distances are larger than the EPM distances 
by ~ 25% on average. The former suggests a value of 
Ho = 55 ± 5. The main uncertainty of this result comes 
from the assumption that £ = 1, where £ is the ratio be- 
tween the total explosion energy and the energy liberated 
by the neutron-proton recombination into 56 Ni (cf. Eq. 10). 
To reduce the P-T distances to the level of the EPM dis- 
tances, which correspond to Ho = 70, an average value of 
£ = 1.9 is required. The consequence that about half of the 
total energy E comes from other sources than the neutron- 
proton recombination into 56 Ni seems rather extreme. In 
fact it is not supported by two SNe IIP (1987A and 1999gi) 
with independent distance information, which suggest that 
£ is of order of unity. Moreover, very recently Leonard et al. 
(2003) have obtained a Cepheid distance of 11.7 ± 1 Mpc to 
NGC 1637 - the host galaxy of SN 1999em, which is by a 
factor of 1.4 larger than the EPM distance (Table 2). Our 
Dp-t distance of 11.1 Mpc to SN 1999em (Table 3) is in a 
good agreement with this result. If it happens that the same 
factor is applicable also to the EPM distances to SNe 1986L 
and 1990E, there will be no need to resort to large ^-values, 
such as £ ~ 2 (section 4), to remove the discrepancy between 
Dp-T and Depm for these SNe. 

In conclusion we emphasize the necessity of construct- 
ing a new grid of hydrodynamic SNIIP models based on 
current evolutionary presupernova models and taking into 
account 56 Ni as an additional parameter in a consistent way. 
Such a grid would allow to create more precise analytic ap- 
proximations for a number of correlations between the phys- 
ical parameters of SNIIP and their observable properties. 

The 'plateau-tail' method of distance determination 
needs, of course, further critical analysis requiring a close 
collaboration between astronomers observing supernovae 
and theorists modelling their explosions. If the proposed E- 
Mm correlation is confirmed it promises to become a tool 
to explore the mechanism of SN II with the aid of optical 
and spectroscopical observations. 



5 CONCLUSIONS 

Model calculation by LN83 and LN85 of SNe IIP, leading 
to Eqs. (l)-(3), are combined with available EPM distances 
and velocity distances (Ho = 60) to derive the explosion en- 
ergy E, the ejected mass M., and the presupernova radius R 
of 14 SNe IIP. Only the apparent, absorption-corrected mag- 
nitude V and the expansion velocity « p h at the mid-point of 
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